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Abstract: We revisit a particular ghost-free bimetric model which is related to both partial 
masslessness (PM) and conformal gravity. Linearly, the model propagates six instead of seven 
degrees of freedom not only around de Sitter but also around flat spacetime. Nonlinearly, 
the equations of motion can be recast in the form of expansions in powers of curvatures, 
and exhibit a remarkable amount of structure. In this form, the equations are shown to 
be invariant under scalar gauge transformations, at least up to six orders in derivatives, 
the lowest order term being a Weyl scaling of the metrics. The terms at two-derivative 
order reproduce the usual PM gauge transformations on de Sitter backgrounds. At the four- 
derivative order, a potential obstruction that could destroy the symmetry is shown to vanish. 
This in turn guarantees the gauge invariance to at least six-orders in derivatives. This is 
equivalent to adding up to 10-derivative corrections to conformal gravity. More generally, we 
outline a procedure for constructing the gauge transformations order by order as an expansion 
in derivatives and comment on the validity and limitations of the procedure. We also discuss 
recent arguments against the existence of a PM gauge symmetry in bimetric theory and show 
that, at least in their present form, they are evaded by the model considered here. Finally, 
we argue that a bimetric approach to PM theory is more promising than one based on the 
existence of a fundamental PM field. 

Keywords: modified gravity, bimetric gravity, Weyl invariance, higher spin fields 



Contents 


1 Motivation and summary of results 1 

2 Review of the derivative expansion in bimetric theory 4 

2.1 The bimetric action 4 

2.2 Perturbative expansion of bimetric equations 5 

2.3 A model with a possible gauge symmetry 7 

3 Linear analysis: PM theory, conformal gravity, and the bimetric model 8 

3.1 Partial masslessness in the Fierz-Pauli framework 8 

3.2 Linearised conformal gravity and PM in bimetric framework 9 

3.3 Linear theory in flat space 12 

4 A systematic study of the extended Weyl symmetry 13 

4.1 General form of symmetry transformations 13 

4.2 Structure of the derivative expansion of the bimetric equations 15 

4.3 Perturbative bootstrap construction of gauge symmetry 16 

4.4 Relation to the partially massless gauge symmetry 20 

4.5 Limitations of the perturbative expansion 21 

5 Discussion 22 

5.1 Summary and discussion of results 22 

5.2 Discussion of the counter-arguments 24 

A Details of bimetric theory and its curvature expansion 27 

A.l Structure of action and equations of motion 27 

A. 2 Obtaining the curvature expansions 28 

B Four-derivative terms in gauge transformations 29 

B. l Higher-curvature equations 29 

B.2 Vanishing of the four-derivative contributions 30 


1 Motivation and snmmary of results 

The two major unresolved issues in General Relativity (GR) are the non-renormalisability of 
the theory and the cosmological constant problem. One expects that a theory of spin-2 fields 
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with more symmetry than GR will be better behaved in these respects. One such theory is 
conformal gravity [1] defined by the action, 

ScG = j ~ J- \R^) (1.1) 

in terms of the Weyl (W) or Ricci (R) tensors; the two formulations differ only by the topo¬ 
logical Euler density. In addition to diffeomorphism invariance this action is invariant under a 
Weyl scaling of the metric. This is a 4-derivative theory and propagates 6 modes, consisting 
of a massless spin-2 field and 4 ghost modes [2-7]. Despite having an extra symmetry as 
compared to GR, the presence of the ghost instability makes the theory less attractive.^ 

Another theory with a novel gauge symmetry is the linear partially massless (PM) theory 
of a massive spin-2 field in an Einstein-de Sitter (EdS) background [11-13]. Due to the 
symmetry, the spin-2 field propagates 4 (instead of the usual 5) polarisations. To this, one 
can add a linear massless graviton (regarding the PM field itself as the graviton requires giving 
up general covariance). The spectrum is now similar to that of conformal gravity linearised 
around EdS backgrounds [7]. However, while ghost-free, this theory is non-interacting and 
exists only around special (EdS) backgrounds. The question is if there exists a unitary theory 
of interacting spin-2 fields with a PM-like gauge symmetry.^ The first investigations involved 
explicit construction of cubic interaction vertices [14-16], but problems are encountered at 
the quartic level. 

On the other hand, a nonlinear setup with the right general features is the ghost-free 
bimetric theory [17]. There exists a unique bimetric model that, when linearised around EdS 
backgrounds, contains a massless and a partially massless spin-2 field, exhibits linear PM 
gauge symmetry and also up to realises the global part of the PM transformations nonlin- 
early [18, 19]. But this is far from establishing the full gauge invariance of the nonlinear 
model. It turns out that the model also exhibits features of conformal gravity. In a derivative 
expansion, it reproduces the conformal gravity equation of motion at the lowest order [20] 
and, hence, is invariant under local Weyl scalings to this order. In this paper we explore 
this property beyond the lowest order and discuss its relevance to the presence or absence of 
PM symmetry. We show that the Weyl scalings can be extended at least up to six-derivative 
terms to maintain the symmetry of the equations. When restricted to the linear theory in 
EdS backgrounds, one recovers the standard PM transformations. The bimetric model pro¬ 
vides a unified description of linear PM theory and conformal gravity up to higher derivative 
corrections. More work is needed to find if this structure can be extended to all orders. 

Another potential setup for PM symmetry is nonlinear massive gravity [21], which is 
ghost-free [22] and can also be formulated around any background, including EdS spacetimes 
[23-25]. A specific massive gravity model, with a fixed dS reference metric, was identified in 

^Nonlinear infinite derivative generalisations geared towards avoiding the ghost problem have been consid¬ 
ered in e.g. [8-10], but the Weyl symmetry is then generically destroyed. 

^To be of relevance to the renormalisation and the cosmological constant problems, the PM symmetry must 
also affect the physical gravitational field. 
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[26] as exhibiting PM symmetry in the decoupling limit. However, it was soon argued that 
this action lacks the gauge symmetry at the nonlinear level [27-29].^ 

There have also been arguments against PM symmetry in the nonlinear bimetric model. 
Since massive gravity models can be regarded as a limit of bimetric models around given 
solutions, it has been argued, for instance in [27, 29], that ruling out PM symmetry in massive 
gravity, implies the same for the bimetric model.^ However, [31] argued that if the bimetric 
model indeed had a PM symmetry, it would most likely be destroyed in the massive gravity 
limit. These results are supported by the hndings in the present work. Other arguments 
are discussed in the last section. Recent work on the subject of partial masslessness includes 
[32-35]. 

Summary of results: This paper does not directly deal with PM symmetry and the pos¬ 
sibility of its realisation in a nonlinear theory. Rather, we explicitly investigate a particular 
bimetric model in a derivative expansion. On eliminating one of the metrics between the 
two equations, one obtains the 4-derivative conformal gravity equation of motion which is 
invariant under Weyl scaling of the metrics [20]. Here we investigate the possibility of ex¬ 
tending the zero-derivative gauge symmetry to higher orders in derivatives. The structure 
of the equations enables us to find a prescription for constructing the gauge transformations 
order by order. 


+ • • ■ , ( 1 - 2 ) 

which leave the equations of motion invariant on-shell up to sixth order in derivatives. Here, 
^{o)9fj.i' = 4>9fj.u and A(o)//ii/ = are Weyl scalings of the metrics and each A( 2 n) contains 
terms with 2n derivatives. The prescription insures that terms with 2n = (4m-|-2) derivatives 
(for integer m) in the transformation can always be constructed. But the existence of the 
2n = 4m-derivative terms requires the model to satisfy certain conditions. It is explicitly 
demonstrated that these conditions are satisfied for the four-derivative term and thus we are 
able to show the invariance of the equations of motion up to six orders in derivatives. This is 
equivalent to considering up to 10-derivative corrections to the 4-derivative Bach equation. 

On Einstein-de Sitter backgrounds, the 2-derivative terms in the transformations of the 
nonlinear metrics reduce to the well-known gauge transformations for a partially massless 
spin-2 perturbation. This suggests that (1.2) may be considered as potential extensions of 
the linear PM gauge symmetry. But wether such a gauge symmetry really exists, remains 
to be seen. In this setup, the metrics are not restricted to de Sitter backgrounds and hence 
we can also study the behaviour of the model around other backgrounds. Unlike its massive 

®The analysis of constraints in [30] rules out massive gravity with 4 polarisations for a Minkowski reference 
metric which is consistent with the linear PM theory. 

"^The massive gravity limit of the bimetric model in [18] does not seem to reproduce the massive gravity 
PM candidate identified in [26], but is consistent with the parameters values found in the subsequent work 
[28]. 
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gravity counterpart, the bimetric model also has flat background solutions which enable us to 
extend the notion of linear partial masslessness to Minkowski space. One finds that around 
flat backgrounds no decomposition into spin -2 mass eigenstates exists, instead, the spectrum 
coincides with the conformal gravity spectrum around flat spacetime, i.e., two massless tensors 
and one massless vector. The results are discussed in more detail in the text and in section 5. 

The paper is organised as follows. In section 2 we review details of ghost-free bimetric 
theory, the perturbative expansion of its equations of motion, and the emergence of the Bach 
equation in a particular bimetric model. Section 3 considers this model at the quadratic 
level and the relation to both linear PM theory and linear conformal gravity, including the 
flat background case. In section 4 we study the perturbative expansions of equations in this 
model and outline the procedure for constructing the higher-derivative terms in the gauge 
transformations of the metrics. We establish the invariance to sixth order in derivatives, and 
discuss the relation to PM transformations, as well as the limitations of the procedure. Our 
results are discussed in section 5, where we also comment on various counter arguments and 
no-go results in the recent literature. Some technical details are provided in the appendices. 


2 Review of the derivative expansion in bimetric theory 

Here we briefly review the structure of ghost-free bimetric theory and outline how to derive 
the higher-curvature expansions of its equations of motion. The main results are summarised 
in this sections; some more details can be found in appendix A. 


2.1 The bimetric action 


The ghost-free action for two spin-2 fields and f^i, with non-derivative interactions is [17], 


S{9J) 





( 2 . 1 ) 


where, nig and mj are the two Planck masses and m is an additional mass scale. The potential 
is given in terms of a square-root matrix S = yg~^ as, 

4 4 

= E (2.2) 

n=0 71=0 

where /?„ are five interaction parameters. The en{S) are the elementary symmetric polyno¬ 
mials of the eigenvalues of 5; for their definitions see appendix A.l.^ 

The fact that both and in (2.1) have Einstein-Hilbert kinetic terms is related to 
the invariance of the potential under the interchanges 5 -H- / and (5n This invariance 

directly follows from the identity [17], 

4 4 

VdY^^neniS) = A/7E^4-ne„(5-'). (2.3) 

71=0 71=0 

®For a non-dynamical = 77 ^,/ and a restricted set of /3„, V{S) becomes the massive gravity potential first 
proposed in [21, 36] and proven to be free of the Boulware-Deser ghost at the nonlinear level in [22, 24, 25]. 
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A Hamiltonian analysis shows that this theory propagates 7 modes [17, 37] and no Boulware- 
Deser ghost [38, 39]. Around backgrounds of the type these modes combine 

into massless and massive spin-2 fluctuations with, respectively, 2 and 5 polarisations [40]. 
In general, a bimetric theory can be interpreted as describing a gravitational metric 
(with standard matter couplings) in the presence of an extra spin-2 field f^y. Then, for 
rrif « rUg = Mp, the physical metric gg_y is a mostly massless field (in contrast to the 
massive gravity limit m/ —>■ oo) [40, 41]. 

2.2 Perturbative expansion of bimetric eqnations 

The equations of motion obtained on varying the action (2.1) with respect to g^y and fg_y are 
of the form. 


g-eom : ^ (i?^, - + ^(5) = 0, (2.4a) 

f-eom : ^ + ^{8-^) = 0. (2.4b) 

In these expressions, Rg_y and R^y are the curvatures of g^y and fg_y, respectively, and Vg_y 
and Vpy are the corresponding interaction contributions, explicitly given in appendix A.l. 
Moreover, S = \Jg~^f. In the g'-equation, the hrst index is raised by g^'^, and in the /- 
equation by To simplify the expressions, we use the notation, 


9 m 

h = — 

nig 


mf 
a = —- 
m„ 


(2.5) 


The ^-equation depends on fg_y (through S) but does not contain derivatives of fg_y. 
Hence, in principle, it can be algebraically solved for f^y in terms of R^y{g)/g?. For example, 
when /3o and /3i are the only non-vanishing parameters, the potential is F = /3o -|- /3iTr(S') 
and its variation gives V^{S) = Pod^y -t- Pi{S^y — 6^yTT{S)). In this case, (2.4a) yields the 
exact expression, 


OM _ I 

I' — 3^1 ^ 






pu •) 


where Ppy is the (scaled) Schouten tensor of g^y, 


( 2 . 6 ) 


R^u — Rfiu g 9^vR ■ 


(2.7) 


Since fg_y = gg_p{S'^Yy, one immediately obtains. 


jtiu — 


Wo 




^T^R/iu T 




Rupgp^Ra 




( 2 . 8 ) 


This is simply a rewriting of the (^-equation, but does not involve solving it as a differential 
equation. Using this in the /-equation (2.4b), gives a higher derivative equation for g^y alone. 

For generic /3n parameters in (2.4a), it is not easy to obtain an exact solution for S, and 
hence for / = gS'^. But, for small Rpy/g?, it is always possible to hnd a perturbative solution 


- 5 ~ 






[20]. At lowest order, i.e., neglecting this gives = a?g^u + • • • j with the constant 

o determined by the polynomial equation = 0. Curvature corrections to this can be 

systematically computed and one arrives at an expression for of the form [20], 

f^iu = o?g^u-\——Z — PPfiu] + O , (2.9) 

where indices on the right-hand side are contracted with g^^. The coefficients a, b and Cn are 
given in terms of bimetric parameters in appendix A.2. Equation (2.9) is a rewriting of the 
^-equation (2.4a) and is satisfied by the compatible solutions of the latter perturbatively,® 
as long as curvatures are small compared to the mass scale ^ = w?/nig. In principle, the 
coefficients in (2.9) can be determined to arbitrary order, and terms with 2n derivatives are 
suppressed by 

Alternatively, we can obtain a solution for g^^i, from the /-equation (2.4b). It has a very 
similar form. 




f j_ p j_ Lp^ j_ z. 

/J, p p 


h2{P)Uu - PP^u 




( 2 . 10 ) 


where Pgu = Pfj,u{f) is the Schouten tensor for Indices are contracted with and the 
suppressing mass scale is now p, = g,/a. 

For generic models, the two expansions (2.9) and (2.10) are not simultaneously valid. 
For example, (2.9) at the lowest order reads = a'^g^u + • • • j whereas (2.10) is of the form 
5 / 21 / = o^// 2 i/ + • • • • Obviously, a necessary condition for the validity of both expansions is 
that (which can be satisfied by fixing one of the (3n). The model we consider in this 

paper satisfies this property. 

It is possible to use the expression (2.9) for to eliminate it from the /-equation (2.4b). 
This yields the following higher derivative equation for g^i,, 


3^oo5/ti/ 4 


I xn p I 3120 R 

i- -pr ngu -I- 






(si + 2S2)P/Ppu - 2s2PPpu - fgpu {PpaPP^ " P^) 


X22 


3PP„„ - 2P/Pp„ - la„„(P^ - P“'’p„a) 




= 0 , 


( 2 . 11 ) 


The coefficients Xmn and are given in appendix A.2, and Qpv = Rpv — ^g^uR is the Einstein 
tensor of gpu- We have collected some of the four-derivative terms into the Bach tensor [1], 

Bpu = -^^Ppu -^p^.P^p + + VpV.P^p - 2P/Pp, + \gp,PP^Ppa . (2.12) 

®Let us comment on the generality of such expansions. If g^i, is an Einstein metric, then 
and (2.9) implies = Pgij.v, for some constant P. On the other hand, in generic bimetric models when one 
metric is Einstein, the other one is Einstein too, but the two metrics are not necessarily proportional to each 
other. Obviously, such non-proportional Einstein solutions are not captured by (2.9). However, in classes of 
bimetric models, the /3i-models (where only one pi out of Pi,P 2 , Pa is non-zero), when either metric is Einstein, 
one necessarily has /an/ = Pg^iv, as implied by (2.9) [31]. In this paper we will work with a /32-model. 


- 6 ~ 









Note that the Bach equation, = 0, is the equation of motion for conformal gravity. 

Equation (2.11) can be rewritten as an Einstein equation for with higher derivative 
corrections. The highest number of derivatives on the same is four; these appear in 
and in some higher order terms which all arise from expanding the Einstein tensor of 
All other terms, including the corrections, contain a maximum of two derivatives on g^i^, 
but in higher powers. This is consistent with the fact that we need to specify four initial 
conditions in the original bimetric equations. 

To summarise, we have re-expressed equations (2.4a) and (2.4b) as (2.9) and (2.11). 
The perturbative equivalence between these two sets of equations holds algebraically and the 
solutions of the earlier set satisfy (2.11) perturbatively (subject to the comment in footnote 
6). For the purpose of this paper, it is not necessary that these two sets are also equivalent 
as differential equations. In particular, studying the symmetry properties of the equations 
involves only algebraic manipulations and, for such purposes, the two sets of equations can 
be treated perturbatively equivalent. 

2.3 A model with a possible gauge symmetry 

Let us consider a particular bimetric model that leads to an equation (2.11) for gfj,i, with, 

a^oo = = 0, (2.13) 

so that the equation starts at fourth order in derivatives. Using the expressions for the Xmm 
it can easily be shown that these conditions uniquely fix the bimetric interaction parameters 
to the following values, 

/3i = /33 = 0, a^/3o = 3a2/32 = (3^ . (2.14) 

It turns out that this choice of parameters also sets X 21 = X 22 = 0 and thus specifies a bimetric 
theory whose equations of motion imply [20], 

+ (2-15) 

At the lowest order, this is the equation of motion for conformal gravity, = 0, which is 
invariant under Weyl scalings g^y —)• since B^j^y —>• e~^^^'>B^y. Hence, at the lowest 

order in a curvature expansion, the equations of the bimetric model (2.14) share the Weyl 
symmetry of conformal gravity. The corresponding transformation of ffj_y can be obtained 
from (2.9). The question is if this is just an accidental symmetry of the bimetric equations at 
the four-derivative level, or if the Weyl scaling could be corrected by adding higher-derivative 
terms to maintain the symmetry at higher orders in the curvature expansion, thereby indi¬ 
cating a gauge symmetry of the full bimetric equations. In this paper we investigate this 
question systematically and, in section 4, show that the corrections to Weyl scaling can be 
computed at least up to six-derivative terms, equivalent to adding up to 10-derivative cor¬ 
rections to the Bach equation. Before that, in the next section, we consider the relation to 
partial masslessness and conformal gravity at the quadratic level. 
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3 Linear analysis: PM theory, conformal gravity, and the bimetric model 

The bimetric model specified by (2.14) is precisely the one identified in [18] based on an 
analysis of partially massless (PM) gauge symmetry in the linearised theory. It propagates 
six modes (instead of the generic seven) around Einstein backgrounds. The question is if it 
can lead to a better understanding of partial masslessness at the nonlinear level. 

Traditionally, PM symmetry is studied in the context of the linear Fierz-Pauli (FP) 
theory. Attempts to find nonlinear generalisations are also often influenced by the FP setup, 
for example, in looking for a fundamental PM held, or in modelling the PM transformations of 
the nonlinear helds after the linear theory. On the other hand, it is also known that conformal 
gravity exhibits a spectrum similar to the linear PM theory, except that now, the underlying 
gauge symmetry is Weyl invariance instead of PM symmetry, and the theory has a spin-2 
ghost. After a brief review of these issues, in this section we consider the bimetric model at 
the quadratic level and show that it provides a unihed description of linear PM theory as well 
as linearised conformal gravity. We argue that the bimetric setup provides a more powerful 
framework for finding a nonlinear generalisation of PM theory. 

3.1 Partial masslessness in the Fierz-Pauli framework 

The Fierz-Pauli equation for a massive spin-2 field in a de Sitter metric is, 

- A{SM^, - + im|p((iM^, - g^,5M) = 0, (3.1) 

where the linearised Einstein operator is given by, 

- V2 + g^^VPV^) 6Mp^ . (3.2) 

It is well known that when the Higuchi bound is satisfied [42] 

TTipp = ^A, (3.3) 

the FP equation becomes invariant under gauge transformations [12], 

= (V^V, + f 5 m.) e(x) • (3.4) 

Consequently, the linear spin-2 field SM^^, has four propagating modes, instead of the usual 
five for a massive field. This is related to the fact that the de Sitter group also has a 
four-component “partially massless” representation, besides the usual hve-component massive 
representation [11-13]. To take gravity into account, we may add a massless spin-2 field 50^^ 
with two propagating modes, to the linear PM theory. 

An interesting question, investigated by many authors [18-20, 26-29, 32, 33, 35, 43], is if 
the linear PM theory in de Sitter background can be extended to arbitrary backgrounds and if 
it can be generalised to a nonlinear theory in a background independent way. However, some 
features of equation (3.1) may be taken as indications that the FP setup is not an adequate 
starting point for such generalisations: 


1. In the FP setup, there is no analogue of PM theory in flat spacetime, indicating a 

preferred choice of background. For A = 0, (3.1) and (3.3) describe a massless spin-2 
field with two propagating modes. In this case, cannot be a PM field since, unlike 

the de Sitter group, the Poincare group has no four-component spin-2 representation. 

2. In a theory with general covariance, cannot be a fluctuation of some background 

independent tensor field This is because (3.1) is not invariant under infinitesimal 

reparameterisations of Then, must the fluctuation of some background 

dependent (and hence, not fundamental) tensor field in such a way that, precisely 
in de Sitter spacetimes, it becomes reparameterisation invariant.^ Below we show that 
these issues can be naturally addressed in a bimetric setup. 

3.2 Linearised conformal gravity and PM in bimetric framework 

Conformal gravity (CG) is defined by the action f d'^x y/g W^, where IF^j/po- is the Weyl 
curvature tensor. The action is invariant under Weyl scalings of the metric g —>■ e'^g. The 
equation of motion is the Bach equation = 0. Due to Weyl invariance, this equation 
propagates six modes of which two are healthy and four are ghosts (or vice versa, depending 
on the overall sign of the action). As observed in [7], in a de Sitter background, the linear 
spectrum is similar to the spectrum of the PM theory with an extra massless spin-2 field, 
except that one of the fields is now a ghost, and the PM symmetry is replaced by Weyl 
invariance. This may be taken as a hint of a connection between CG and PM theories. An 
attempt to identify a nonlinear PM field in CG was made in [43] but no such field was found. 
Note that unlike the PM theory in the FP framework, conformal gravity admits flat space 
as a background around which it propagates six modes; two massless spin-2 fields and one 
massless vector [6]. 

Now, we consider the bimetric model (2.14) at the quadratic level and show that it 
provides a unihed description of both linear PM theory and linearised conformal gravity. The 
Einstein-de Sitter solutions in this model are of the type = c^g^y where the equations leave 
the constant (? arbitrary. The arbitrariness is unique to this model and is a consequence of 
PM symmetry (see footnote 8 below). In these backgrounds, the fluctuations 5g^j,y = g^y—g^y 
and Sffj,y = ff^y — f^y combine into massless and massive spin-2 modes [40], 

dG I^y dg^y -|- (y. df^y , dM^y dffj^y C dg^y . (3.5) 

The cosmological constant and Fierz-Pauli mass in this model satisfy the Higuchi bound and 
are given by (in the background metric g^y and using the notation in (2.5)), 

Kg = |m|p = 3/32 + c^) • (3.6) 

^This leaves out the possibility that is a nonlinear field with some nonlinear completion of (3.1). 
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The linearised bimetric action, diagonalised into the above mass eigenstates is, 


C _ 
^lin — 


J d^x - ^{6G>^’'SGf,, - ^dG^) (3.7) 


where £ is the linear Einstein operator defined in (3.2). This action is invariant under lin¬ 
earised diffeomorphisms of 5Gpiy as well as the linear PM gauge transformation, 

A{5Mp,) = A{Vpd, + ^gp,)C{x), (3.8) 


with gauge parameter ^(x) and for any A. In principle, this may also be accompanied by a 
(restricted) coordinate transformation of SG^y, 


A{5Gpy) = BVpdyi{x). 


(3.9) 


Using (3.5) one can easily work out the corresponding transformations for the original vari¬ 
ables A{5gpy) and A{6fpy).^ Either A or B can be absorbed in a rescaling of as long as 
it remains non-singular. Beyond that, the arbitrariness in A and B cannot be fixed in the 
quadratic theory. In section 4.4 we show one can set B = 1 and nonlinear considerations in 
bimetric theory determine. 


A = 


1 

2c^ 


(1 — a^c^) 


(3.10) 


Owing to the above gauge symmetries, SM^y is a partially massless field and SG^y is a 
massless field, with four and two propagating modes, respectively. Eor positive values of c^, 
non of the modes is a ghost and (3.7) is a healthy action. 

As pointed out in [7], conformal gravity in a de Sitter background has the same structure 
of modes as above, except that the PM field is a ghost. The reason for this similarity can be 
easily understood at the quadratic level. For negative values of in the action (3.7) the field 
6Mpy becomes a ghost because the sign of its kinetic term changes. In this case the linearised 
bimetric action can be related to the linearised conformal gravity action. To see this, consider 
0 > > —a^, and rescale the fluctuations by real constants (with |c| = V—c^), 


<5M, 




9 arug 




followed by an additional field redefinition. 



(3.11) 


6G'py = 5Gpy - ^6Mpy . (3.12) 

®As shown in [18], for constant these transformations can be integrated to hnite ones only if is 
undetermined. This, in turn, uniquely led to the parameters (2.14). Thus, changes in are associated 
with constant PM “gauge” transformations restricted to the proportional backgrounds. For more general 
cosmological backgrounds, the equations leave a time-dependent function undetermined [18]. 
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Only for (? < 0, this replaces the kinetic term of 5M^i, by a kinetic mixing term, 

> ( - ^ {5G'^^^5G'^, - \6G'^) ) 

- 26M^^£^'^P^6G'p^ + Ag {6MP’^6G'g^ - ^6M6G') + 5M^^''5Mg^ - 6M^ 


S' = ^ f 

'^Im 1+a^c^ J 


d^x 


(3.13) 


This action is the linearised form of the auxiliary-field formulation of conformal gravity [3], 


<? — 


TW / 


^ (i? - 2A,) - ( Gg, + AaG'„„ ] + - M 


9^ fiu 




(3.14) 


where R and Qgu are the scalar curvature and Einstein tensor of G',jy. The solution to the 


(lU • 


A'lgv equation of motion, 


= -^G'g, + \ [Rg, - \G'g,R) , 

when plugged back into (3.14), gives the conformal gravity action for 


(3.15) 


5c = -id+k^ / {R^''R^.u - \R^) ■ (3.16) 

It is then obvious that integrating out in (3.13) leads to the linearised form of this 

action. We emphasise that the manipulations (3.11) and (3.12), remove the kinetic 

term only if c2<0. Thus, considering the bimetric model in a regime with spin-2 ghosts is 
crucial for making the connection to conformal gravity at the level of the action. 

We conclude that bimetric theory, around its proportional backgrounds, gives a unified 
description of linearised conformal gravity and linear PM theory. At the quadratic level, both 
these theories are different phases of the same bimetric model. It should be emphasised that 
although the undetermined modulus parameter is equivalent to a gauge parameter, it is 
not possible to start with > 0 and reach < 0 by continuous gauge transformations. This 
is because such transformations would necessarily have to cross = 0, for which one of the 
metrics becomes singular. Hence, the ghost-free linearised PM theory does not lie on the 
same gauge orbit as the linearised conformal gravity. 

As a final remark, the quadratic bimetric and CG actions are equivalent only for (? < 0, 
while for (? > 0, the actions are not equivalent. However, for both > 0 and (? < 0, the 
linearised bimetric equations lead to the linearised Bach equation. For any c^, the linearised 
bimetric equations are, 

^ {Sgpu - - c^5gg^ - {6f - c^Sg)ggu^ = 0, (3.17) 

^ gpy) + 3(i+t^e^) - c^^9pu - (5/ - c^5g)ggu'^ = 0 . (3.18) 

Solving the first of these equations for 6fgi, and plugging the result into the second always 
gives the linearised Bach equation for 5gfj_,y, 

^6g^0 = O. (3.19) 

0gal3 
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3.3 Linear theory in flat space 


The discussion in the previous subsection is valid only for Ag ^ 0. However, unlike the PM 
theory in the FP formulation (3.1), the bimetric model also has a flat space solution with 
six propagating modes around it, without a conflict with the Poincare group representations. 
From the expression for Ag (3.6) it is evident that the flat space solution corresponds to a 
choice of 


= -a-^ . (3.20) 

As mentioned above, this value cannot be reached by continuous transformations from > 0.® 
With the choice (3.20), the expressions of the massless and massive modes in (3.5) coincide 
or, equivalently, the expressions for and 5fgy in terms of the mass eigenstates become 
singular. It is therefore not possible to diagonalise the equations in terms of spin-2 mass 
eigenstates. These are no longer the correct variables to work with, and the earlier discussion 
of the relation to conformal gravity, formulated in terms of 5Ggi, and 6Mgy, is not valid. This 
also resolves the conflict with the absence of spin-2 PM representations of the Poincare group. 
The equations are, however, well-defined in terms of the original and 6fgy and become 
most transparent when expressed in terms of dg^y and 6M^y = 6Ggy = 5g^y + oP'hf^y, 

= 0 , (3.21a) 

£'gy5gpa + V'/32 [^Mgy - ^gybAlP) = 0 , (3.21b) 

where E'^^y is the linearised Einstein operator (3.2) in flat space. Despite the form of the first 
equation, 5Mpy is not a massless spin-2 field because the entire system is not invariant under 
linearised coordinate transformations of 5Mpy. Coordinate transformations of hg^y, 

^ ^9pu T dg^y dy^fi , (3.22) 

are symmetries of (3.21) but dg^y does not satisfy the equation of a massless spin-2 field. 
The system of perturbation equations is also invariant under a scalar gauge symmetry which 
transforms both dg^y and 5Mgy, 

A6ggy = gpy , A6Mgy = dgdyC - (3.23) 

Moreover, taking the divergence of (3.21b) shows that 6Mpy is transverse. Therefore, the 
number of propagating degrees of freedom around flat background is 

20 components — 2-4 coordinate transformations 

— 2-1 scalar gauge transformations — 4 constraints = 6. (3.24) 

®Flat space solutions do not exist in the corresponding massive gravity model of [28]. In this case, the 
cosmological constant of the proportional backgrounds is Ag = where again c is undetermined. 

Implementing Ag = 0 by setting = 0 leads to a singular metric (since A,/ is fixed). More generally, in the 
massive gravity limit, the CG-related phase of the bimetric model as well as many other solutions disappear 
[31], although more pathological solutions for < 0 could survive. 
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This is the same number of modes as found around a de Sitter background, where the fluc¬ 
tuations are diagonalisable into a massless and a partially massless field. Around flat space, 
however, no interpretation in terms of mass eigenstates with spin-2 alone exist. In fact, the 
spectrum is the same as that of conformal gravity which consists of two massless spin-2 and 
one massless spin-1 fields [2, 4-6]. Depending on the overall sign of the action, four or two 
out of the six propagating modes are ghosts. 


4 A systematic study of the extended Weyl symmetry 

Motivated by the appearance of the Weyl invariant Bach equation at lowest order in the 
derivative expansion of bimetric equations (2.15), we now investigate the possibility of ex¬ 
tending this invariance to higher-derivative terms in the expansion. 


4.1 General form of symmetry transformations 

Let us consider the bimetric model (2.14) in more detail. The action is given by (2.1) with 
the potential. 


y = 3^2 (a-2 + 1{[S^] - [5]2) + ^ 

where [S'] = Tr 5 and S = -s/g~^f. The g and / equations of motion are, 

+ (S^ - = 0, (4.2a) 

^i?^,(/) - 3^ + ^(S-2 - [S-^]S-^r, = 0, (4.2b) 

in which we have raised the first index with g^'^ and respectively, and used the notation 
in (2.5). Note that in terms of rescaled variables, 

9fj,u ~ ^ 9iJ.y ■> — oif ) ( 4 :- 3 ) 

the equations (4.2) take a very symmetric form, 

aR^,{g') - 350 + = 0, (4.4a) 

aR^W) - 350 + {S'-^ - [S'-^]S'-^)0 = 0, (4.4b) 

in which we have used the further notation (see (2.5)) 


These transform into each other under the interchange g'^^, o keeping all parameters 
fixed. This property will play an important role in the considerations below. 

Our aim is to check if we could find nontrivial, field-dependent transformations g' —>• 
g' + Ag'{g',f') and f ^ f + Af'{g',f') (involving also a scalar gauge parameter) that 
keep the above equations invariant. Since dealing with the variation 5S of a square-root 
matrix is cumbersome, we do not directly study the variations of (4.4), but rather attempt 
to construct the symmetry perturbatively. First, we deduce some general properties of such 
transformations. 
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Interchange symmetry: The ^''-equation (4.4a) can be rendered invariant for any Ag'{g', f) 
by choosing an appropriate compensating transformation Af'(g', f). Such a Af can be com¬ 
puted for any Ag' at least perturbatively, for example, by using the perturbative solution for 
/' in terms of g' and its curvatures in (2.9). To find a symmetry of both equations, one 
must be able to specify Ag' such that the same variations also keep the /'-equation (4.4b) 
invariant. Let us assume that such transformations exist. 

On the other hand, since the /'-equation (4.4b) is obtained from the ^'-equation by simply 
interchanging g' and /', it must also inherit the symmetry of the ^'-equation. In other words, 
the /'-equation must also be invariant under, 

^new5 ~ ^/ l/'-O-g' j ^new/ = Ag . (4.6) 

A similar argument holds for a new symmetry of the gf'-equation. But if we know from 
various limits that the equations can at most admit one symmetry, then the two sets of 
transformations must coincide, Ag' = A^ewg', A/' = Anew/^ that is, 

^ 9 ' = ^/l/'-H-g' ) ■ (4-7) 

This property will become very useful in constructing the symmetry transformations. 

Field dependence convention: The appearance of the square-root matrix, S = \Jg~^f, 
in the above equations makes it difficult to investigate their symmetry directly by treating 
the fields non-perturbatively (although still linearly in the symmetry parameter). Hence here 
we will work with the derivative expansion of the above equations by recasting the ^'-and 
/'-equations as the expansions (2.9) and (2.10) for /' and g', respectively (see below). Then 
the gauge transformations can also be expanded as, 

Ag' = ^ , A/' = ^ A( 2 „)/' , (4.8) 

n=0 n=0 

where the variation A( 2 n) contains a total of 2n derivatives of the fields and the gauge param¬ 
eter. In this approach, the expression Ag'{g',f') can be expanded in different ways since a 
dependence on /' can be re-expressed in terms of g' using (2.9), and vice versa, thereby mix¬ 
ing different orders of curvatures. This ambiguity can be avoided by following the convention 
that Ag' is expressed entirely in terms of g' and its curvatures, such that Ag' = Ag'{g') and, 
similarly, Af = A/'(/'). 

The bootstrap construction: In the previous section we saw that, on perturbatively 
eliminating /' between equations (4.4), one obtains the Bach equation for g', plus corrections 
involving higher powers of Rg_u{g')/9^■ Thus, at the lowest order, A^^Q^g'^^^, = is a sym¬ 
metry. To check if this symmetry extends to higher orders, a straightforward approach would 
be compute the higher derivative corrections to the Bach equation and see if corresponding 
higher order corrections A( 2 n) 9 ' to the symmetry transformation exit. This would show the 
existence of an on-shell symmetry in this perturbative framework. The variation Af could 
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then be computed by plugging lS.g' in the expansion (2.9) and expressing the outcome entirely 
in terms of f' using (2.10). The A.g' and A/' obtained in this way will automatically satisfy 
(4.7). In this approach, to compute, for example, the 4-derivative term in lS.g' one needs to 
know the 8-derivative correction to the Bach equation, which is rather lengthy to compute 
and manipulate. 

Below we follow an alternative bootstrap approach where, starting with A(o)g'^ the trans¬ 
formation is systematically constructed order by order, by using the perturbative expansions 
for /' and g' and imposing the interchange symmetry (4.7). The transformations are then 
guaranteed to leave (4.4) invariant order by order. 

4.2 Structure of the derivative expansion of the bimetric equations 

As discussed in section 2.2, the equation (4.2a) contains algebraically through S^y 
and hence can be solved perturbatively to express in terms of R^y{g)/g?I32- A similar 
statement holds for the f^y equation (4.2b).^® For the model (4.1), the expressions (2.9) for 
ffj,y and (2.10) for g^y obtained in this way become, 

CX) 

fii.u = + 4 X] ’ (4.10a) 

n=2 

oo 

+ 4 ^ [/] ■ (4-10b) 

n=2 

In terms of the rescaled variables (4.3) and (4.5), these take an interchange symmetric form, 

OO 

+ aP,, + 4Y.(l)" T"’[9'I . {'l-lla) 

n=2 

oo 

+ &P,^ + 4E (t)” • (4.11b) 

n=2 

Here ^^^v\g\ denote terms with 2n derivatives, involving n powers of P^y and (n — 1) powers 
of the inverse metric g^’-' to contract the indices. In appendix B.l, we provide their explicit 
form for terms up to fifth order in P^y, i.e. up to n = 5. The same functions appear in 
both equations due to the g' -f-)- f interchange symmetry of the model. Starting from (4.4), 
one can easily write down recurrence relations that generate to any order. Note that 
such expansions are not valid in the massive gravity limit {a —>■ oo) of bimetric theory. 

other words, equation (4.4a) can be viewed as algebraically determining S' in terms of the matrix 
aR'^y{g') — 35Ti or aP'Kg') — after subtracting jR. So, in principle, (4.4a) can be re-expressed as, 

S'^ F {&P>i{g') - , (4.9) 

where F denotes some matrix function. Since P't commutes with (5A, for small enough aP'l, F can be 
expanded in a power series in the usual way. Alternatively, this power series can be obtained directly by 
perturbatively solving (4.4a) for S', which is our approach here. Finally we compute (4.11a) from / = gS^. 
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Equations (4.11) are perturbatively equivalent to the equations of motion (4.4) in the 
algebraic sense (for the purposes of symmetry arguments that involve algebraic manipulations 
alone, it is not necessary that they are also equivalent as differential equations). Also, at the 
lowest orders, the two expansions are compatible.Finally, we note that all 7 ^^”^ possess 
the following important properties: 

• vanishes identically on Einstein backgrounds for which P^i, oc g^,y. 

• Beyond that, 7 ^^”^ vanishes identically on the homogeneous and isotropic solutions 

of [44] for which, in the above model, oc + a^/^i/). 

• Linear perturbations of 7 ^?/”^ around Einstein backgrounds also vanish identically. 


.ST (2n) 
^9 pa 


^9pa — 0 • 

Einstein 


(4.12) 


These statements can straightforwardly be proven by noting that for the above type of back¬ 
grounds, the complete equations of motion imply = I52p^{gpv + cP'fpu)- But this is 
saturated by the first two terms in the above expansions, hence the contributions from higher 
terms must vanish.This implies that there is no contribution from any of the 7 ^^”^ on these 
backgrounds. In the same way one verifies the statement about the linear variations. 

The above properties imply that for backgrounds of the type / = c^g, all 7 ^^”) vanish 
and the expansions in (4.10) or (4.11) become exact. At the lowest order, P^y = 0 and one 
gets / = —g. At the quadratic order, P^^y = {Ag/3)gfj,y leads to / = cPg, with no higher 
derivative contributions (with cP and Ag related by (3.6)). At the end of this section we will 
comment more on the validity regime of these expansions. 


4.3 Perturbative bootstrap construction of gauge symmetry 

Now, employing the criteria collected in section 4.1, we find the condition under which a 
gauge symmetry could exist and describe a procedure for constructing the transformations 
to any order. We explicitly show that a gauge symmetry exists to, at least, six orders in 
derivatives. 

Zero-derivative terms: Guided by the Weyl invariance of the Bach equation, which ap¬ 
peared in the analysis of section 2.2, we take the lowest-order gauge transformation of (4.4) to 
involve an infinitesimal Weyl scaling of g'. Then, the g' equation (4.4a) can always be rendered 

^^The purpose of these expansions here is not to accurately compute perturbative corrections, in which 
case one has to ensure that the two curvature expansions are mutually compatible. The idea rather is that if 
equations (4.4) admit a gauge symmetry, then the expanded equations (4.11) will exhibit a symmetry of the 
form described here, irrespective of the details of the compatibility of the two curvature expansions. 

^^Up to n = 5, this can be verified directly from the explicit expressions provided in Appendix B.l. 
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invariant by a compensating transformation of f which can be easily computed perturba- 
tively using the expression (4.11a) for /' in terms of g'. To lowest order in derivatives, this 
gives, 


Using (4.11b), A(o)/(n/ can be expressed in terms of alone. Thus to zeroth order in 
derivatives, one gets. 


\o)9'^u = 4>9'^^u, = Pf^u- (4-14) 

Obviously, the lowest-order terms in the transformations satisfy the criterion of interchange 
symmetry (4.7) and hence must keep both equations (4.4) invariant to this order. It is trivial 
to see that on ignoring derivative terms the interaction contributions in (4.4) are invariant 
under such transformations.^^ 

Two-derivative terms: Now we describe in detail the construction of the two-derivative 
terms in the transformation. The procedure generalises to higher orders straightforwardly. 

• Step 1: The first step is to determine the structure of all two-derivative terms that are 
generated by the zero-derivative part of the transformations. Using A(o)5^ from (4.14) in 
(4.11a) and retaining terms with up to two derivatives gives the following compensating 
transformation. 


+ ^{2) flu = -pglu - aV^dup + ... . (4.15) 

The dots denote still missing two-derivative terms. The field dependence ambiguity 
discussed in section 4.1 is fixed by re-expressing Af entirely in terms of f. Using the 
expression (4.11b) for g' in terms of ff the above transformation at the two-derivative 
level then becomes (pfl^, — Oi(j)Pfj,u — aV^duf- Since the equations treat g' and /' on 
equal footing, it follows that similar two-derivative terms must also appear in A^ 2 ) 9 '■ 
Hence we take, 

\o)9lu + ^{2)9lu = P9lu + aiPPiiu + a2V^d^(j), (4.16) 

with constants oi, 02 to be determined. Now, using this in (4.11a) to recompute A/))^,, 
and once again expressing its ^'-dependence in terms of /' gives, 

\o)flu + ^{2} flu = Pflu - (oi + d) (j)Pf,u - (02 + d) Vf,du(j). (4.17) 

By construction, these leave (4.11a) invariant on-shell up to second order in derivatives. 
No new two-derivative terms are generated. 

course, the interaction contributions in (4.4) are invariant under g' —>■ Ag' , f' —> Af' , for any invertible 
matrix A, but only the Weyl part of this admits an extension to include derivative terms. 
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• Step 2: To ensure that the same transformations also keep (4.11b) invariant, it is suf¬ 
ficient to impose the g' ■H' f' interchange symmetry (4.7), as argued in section 4.1. 
Comparing (4.16) and (4.17), one sees that this requires oi = 02 = —d/2. Hence 
the transformations that leave both equations in (4.11) or (4.4) invariant at the two- 
derivative level are, 

- f , (4.18a) 

A(o)/;. + A(2)/;, = </>/;, - f . (4.18b) 

Let us briefly comment on some ambiguities at the two-derivative level that are avoided by 
the above construction. One can check explicitly that even for arbitrary ai and 02 , the 
transformations are a symmetry of (4.11) at the two-derivative level (although this is fixed 
by the higher-derivative terms). This is due to the fact that at the two-derivative level, 
a coordinate transformation 6x^ = —a 2 g'^''d^ 4 > induces variations 5g'^^^ = 2a2'V and 
= —2o2V^9i/(/>. Also, to this order, the equations (4.11) are unchanged under field 
redefinitions g' ^ g' + aiP and f ^ f — aiP. Hence, ai and 02 are not genuine parameters 
of the new symmetry. In this case, imposing the exchange symmetry (4.7) not only guarantees 
the invariance of (4.11b), given that (4.11a) is invariant by construction, but it also fixes the 
coordinate and field redefinition ambiguities. 

Extension to higher orders: The construction of the two-derivative terms straightfor¬ 
wardly generalises to higher orders. Two situations arise depending on the number of deriva¬ 
tives (2n -|- 2) for even or odd n. Given transformations Ag' and A/' to 4m-derivative order 
(odd n), they can always be extended to (4m -|- 2)-derivative order. But starting with Ag' 
and Af to (4m — 2)-derivative order (even re), an extension to 4m-derivative order exists 
only under certain conditions. These need to be satisfied in order for a complete perturbative 
expression of the gauge symmetry to exist. Remarkably, the four-derivative terms turn out 
to meet this condition and hence we are able to prove the presence of a gauge symmetry up 
to six orders in derivatives. 

In order to demonstrate the above, let us try to generalise the two-derivative construction 
to terms with arbitrary number of derivatives in the spirit of a proof by induction. For this 
we assume that X]fc=o A( 2 fc)<?^ a-iid Yllt=o^{ 2 k)f' are gauge symmetries of (4.11) at the 2re- 
derivative level and satisfy the interchange symmetry (4.7). This implies that if we use 
^^^gA( 2 fc) 5 ^ in (4.11a) and express the resulting variation for f entirely in terms of f 
using (4.11b), then we get Ylk=o ^{ 2 k)f'■ Moreover, the interchange symmetry insures the 
invariance of (4.11b). 

We are now looking for transformations ^{ 2 k) 9 ' and ^( 2 k)f' that keep the 

equations invariant to next order. The new purely (2re -|- 2)-derivative terms, A( 2 n-i- 2 ) 5 ^ and 
A( 2 „+ 2 )/^ can be constructed following the two steps outlined earlier. 

• Step 1: We start by finding all relevant (2re -|- 2)-derivative terms, i.e., those that are 
generated from lower orders. Using Ylk=o ^{2k)9' ia (4.11a) to compute the variation 
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for /' produces directly a set of (2n + 2)-derivative terms in terms of g'. Furthermore, 
converting all g' into f using (4.11b) produces additional terms that will appear in 
A( 2 n+ 2 )/^-^^ Motivated by the interchange symmetry (4.7), we then introduce all the 
corresponding terms with arbitrary coefficients into A( 2 n+ 2 ) 5 ^ and recompute the form 
of A( 2 n+ 2 )/^ from (4.11a). As a representative of the ( 2 re + 2 )-derivative terms in the 
transformations, let us consider, for instance, a term of the type (for a constant 

A) which is generated in Af on using Ag' at order 2n in (4.11a), and expressing the 
result in terms of g'. This term contains n powers of to contract the indices. The 
corresponding term must be included in A( 2 „_|_ 2 ) 5 ^ with an arbitrary coefficient B, 

A(2n+2)V = i?<^/^;.+' + -- - • (4.19) 

Now when A( 2 n+ 2 )/^ is recomputed from (4.11a) with A(^2n+2)9' included, one gets, 

A(2n+2)/;. = ((-1)"+'S + A) + • • • . (4.20) 

The A-term is generated from lower orders, while the S-term arises directly from (4.19). 
The factor of (—l)"'"''^ arises from the replacement g'^^ = + ... in (we pick 

up a minus sign for each of the n factors of g'^^), along with the overall minus sign 
for computing the transformation of f' from the first term of (4.11a). By construction, 
these transformations leave (4.11a) invariant. 

• Step 2: The second step is to impose the interchange symmetry (4.7) in order to insure 
the invariance of (4.11b). This yields the condition, for any single term, 

B = {-ip^^B-A. (4.21) 

Clearly, this can be solved for B only for even n = 2m. This means that we can impose 
the interchange symmetry guaranteeing the invariance of both (4.11a) and (4.11b) only 
for terms with (4m + 2) derivatives. For terms with odd n = 2m — 1, or 4m derivatives, 
the construction fails unless A = 0. This must hold for all types of terms generated at 
this order. If this condition is met, i.e. if all contributions generated from lower orders 
cancel each other out, then the coefficients B are left arbitrary by the construction at 
this order in the transformation, though they may be fixed by higher orders. 

We conclude that an on-shell gauge symmetry can be constructed provided that, at each 4m- 
derivative level, the contributions generated from lower orders vanish on-shell. (Dis)proving 
the existence of a nonlinear gauge symmetry in the equations (4.11) or (4.4) reduces to 
checking whether this condition is violated or not. 

^^Note that all terms with lower number of derivatives are already accounted for in Y^k^o A{ 2 k)f'■ 
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Four- and six-derivative terms: The above analysis shows that, in this approach, an 
obstrnction first arises at the fonr-derivative level. For a gange symmetry to exist, it is 
necessary that the variation A(o)5' + ^( 2 ) 9 ' used in (4.11a) does not generate any fonr- 
derivative terms in A/'. This is not obvious at first sight. However, the outcome of a rather 
lengthy calculation outlined in appendix B is that such four-derivative terms are indeed not 
generated by terms with fewer derivatives (this can be easily verified for a constant gauge 
parameter (p). Therefore, the four-derivative terms in the transformations remain arbitrary 
at this level. They may need to be fixed when considering the eighth-order contributions, 
which again need to vanish on-shell due to (4.21). 

The above construction also shows that, once a symmetry at the four-derivative level 
is established, invariance at the six-derivative level is guaranteed. Hence we are able to 
perturbatively construct gauge transformations of (4.11) up to sixth order in derivatives, 

^9'fj,u = ‘^{0)9'fj,u + + ^{a)9'^iv + ^{6)9'^iu + • • • , (4.22a) 

+ • • • ■ (4.22b) 

The result for the zero- and two-derivative terms is given in (4.18). Terms with four deriva¬ 
tives in the transformations are completely arbitrary at this stage; any choice will leave the 
equations invariant up to sixth order in derivatives. The six-derivative terms, and 

A(g)/', are again determined by the construction described above (but are expected to de¬ 
pend on the choice of the four-derivative terms). The explicit expressions are tedious, but 
the derivation is straightforward.^^ We emphasise that these transformations are covariant 
and background independent to the extent that the perturbative expansion is valid. 


4.4 Relation to the partially massless gauge symmetry 

In section 3.2 we described how the linear PM theory is easily embedded in a covariant 
bimetric setup. In this section, we have constructed transformations (4.22), at least up to 6 
orders in derivatives, that keep the bimetric equations invariant. In the perturbative setup 
employed, the transformations are simply generalisations of the Weyl scalings of the metrics 
augmented by higher derivative corrections. This construction required no input from the 
PM properties of the quadratic theory. 

It is now easy to verify that the linear PM transformations (3.8), (3.9) are recovered from 
(4.22) at the two-derivative level. This also fixes the arbitrariness of the quadratic theory. In 
terms of the unsealed metrics g and /, the two-derivative variations (4.18) can be rewritten 
as (after adding and subtracting Kg contributions using (3.6)), 


^giMU ~ ^ )4‘9fiu 2 / j ,^/32 

= ^(1 - - 2 ^ 


P — 

fiiy 3 


p _ f 
flu 3 ^ J 


+ /idu 


fl'-'UY I J 


(4.23a) 

(4.23b) 


°To find these transformations using corrections to the Bach equation, we would have had to compute up 
to 10-derivative corrections to the Bach equation. 
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where, higher derivative terms have been suppressed. Let us restrict the nonlinear fields 
and to g^y = g^^y + dg^y and f^y = f^y + 5f^y. The fixed background satisfies ff,y = c^g^^y, 
so the transformation only affects the fluctuations, /S.g^y = A{6g^y) and Af^y = A{6ffj,y). 
For such a background, P^y = P^y = {Ag/3)gfj_y. Since the transformations are already linear 
in the small gauge parameter cj), only the background pieces contribute to the right-hand-sides 
of (4.23). Then, at two-derivative level, the symmetry transformations reduce to, 

’ (4.24a) 

A{5fgy) = -^ (1 - a^c^) ^g^.u - (4.24b) 

For the mass eigenstates (3.5), this implies, 

A{6Mgy) = (V^V, + ^-gg^ </-, (4.25a) 

A(<5G^.) = -^VgVyC^ (4.25b) 

This recovers the standard PM transformation of the quadratic theory in Einstein-de Sitter 
backgrounds, (3.8), along with a coordinate transformation of (3.9) (for a^cf) = —g?‘l320- 
also fixes the ambiguity present in the quadratic theory. Note that AM^y is invariant under 
infinitesimal coordinate transformations of g and /. 

In this approach, it is obvious that the standard form of the PM transformation in terms 
of 6Mgy is very specific to Einstein-de Sitter backgrounds. Away from such backgrounds, the 
fundamental fields are g^^y and f^y, and a fundamental PM field does not exist. Purthermore, 
the nonlinear transformations induce no linear variation of SM^y when = 1. We will 
comment on an implication of this in the next section. 

4.5 Limitations of the perturbative expansion 

The perturbative expansions we have used are valid for small Rgy{g) / (/ 32 /i^), or small aPgy{g'), 
where l// 32 /i^ or d are the small expansion parameters. To connect to the linear PM theory, we 
considered the expansions around Einstein-de Sitter (EdS) solutions for which f' = a^c^g'. 
But these backgrounds have a curvature scale A'g ~ l/d (as implied by (3.6) in terms of 
rescaled variables). Thus, in these cases, 

aPfj.y{g') = (1 + a^c^)gfj,y (4.26) 

is not suppressed by a and is small only when (? is close to In spite of this, the standard 

PM results were recovered for any (? from the first two terms in the expansions. The reason is 
that although now the 4- and higher-derivative terms 'yj^J^\g') are not suppressed by a, they 
nonetheless vanish on EdS backgrounds and, hence, are small in the vicinity of it, for reasons 
having to do with the structure of the equations. The linear variations of 7 ^^”^ around EdS 
solutions also vanish. But, now the expansion in 7 ^^”^ is not perturbative since terms with 
different n may make small, but comparable contributions. Hence the expansions we have 
used are definitely reliable as perturbative expansions when |1 -|- a^c^| < 1, in which case 


- 21 - 






the spectrum of perturbations is not unitary. However, if the presence of a symmetry could 
be established to all orders in this regime, this would also imply the invariance of the full 
equations of motion. For c2>0, the status of the expansion is less obvious although the linear 
results around EdS solutions are easily reproduced. In any case, the possibility of extending 
the symmetry beyond the four-derivative Bach equation clearly points to additional structure 
within the nonlinear equations. 

It is easy to rewrite equations (4.4) such that they are valid perturbative expansions 
around general proportional backgrounds, but this is not useful for finding the transforma¬ 
tions. In (4.4) we can simply rewrite R^y{g') as {R^y{g') — h'gg'^^) + R'gg'^y and expand S' in 
powers of a{R^y{g')-K'gg'g^), or equivalently, in powers of aP^y{g') = a{P^y{g')-{K'g/?>)g'g^). 
In the notation of footnote 10, the equations imply that S' must have the form, 

S' = F + (d^ - = F , (4.27) 

where F{X) is a matrix function of the matrix X and traces of powers of X. Since 
commutes with for small enough d, this can be expanded in a power series in the usual 
way. In practice, this is again achieved by solving the equation perturbatively for S'. The two 
lowest-order terms are just rewritings of the corresponding terms in (4.11), = a^c^g'^^ -|- 

aVfxy{g') -!-•••. At higher orders, the are now functions of V, rather than F, and the 
coefficients are c-dependent. Close to EdS backgrounds, aV^u are small and the expansion 
is perturbative. However, under a variation, aXV^u = aXF^y — (1 -|- a^c^)Ag'g^, note that 
Ag' is not suppressed by a small d. Hence higher terms contribute to the variation of lower 
terms and the expansion is not useful for finding the transformations order by order. 
Alternatively, in (4.4) we can rewrite Fgy{g') as [F^y — ^A'gg'g^y) + ^A'gg'gy. Then, 

S'= F (dn + (1 + «2c2) g>^%^ - 5^) . (4.28) 

Now, taking = 0, the variation of V^y = F^y — ^A'gg'^^y does not mix terms of different 
order. But F^y does not commute with {g'~^g')^u and the expansion is involved. Hence there 
is no easy way of improving over the expansions considered here. 

5 Discussion 

In this section we first discuss our results and then argue that the bimetric model considered 
here avoids the counter arguments presented so far for the absence of a gauge symmetry in 
the bimetric setup, at least in their present form. 

5.1 Summary and discussion of resnlts 

In this paper we investigated a possible gauge symmetry of the equations of motion of a 
particular bimetric model specified by (2.14). The equations contain the square-root matrix 
S = \/9~^ f s-iid the expression for the variation of S in terms of the variations Ag and 
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A/ is complicated. We avoid this complication by working with appropriate perturbative 
expansions of the equation. 

In a small curvature regime, we perturbatively eliminate one of the metrics, say, 
between the two equations to get a higher derivative equation for g^y. At the lowest order, 
instead of the Einstein equation for for this particular model one obtains the 4-derivative 
Bach equation of conformal gravity which is invariant under Weyl scalings of In a 

systematic treatment, we show that the Weyl scalings of g^y and f^y can be corrected at 
least up to 6-derivative terms to maintain the gauge symmetry at higher orders. A brute 
force computation of the same transformations, by demanding the invariance of the higher- 
derivative extension of the Bach equation, would have required calculating up to 10-derivative 
corrections to the Bach equation. Furthermore, we find a sufficiency condition for extending 
the symmetry to higher orders: In our approach, at orders 4m there are obstructions to the 
symmetry which must vanish. Then a symmetry at levels 4m -|- 2 is insured to exist. Once an 
obstruction vanishes at a given order 4m, the transformations at that level remain completely 
arbitrary and will be fixed at order 4(m -|- 1). 

The above results are completely independent of de Sitter backgrounds and linear partial 
masslessness. Nevertheless, it turns out that in Einstein-de Sitter spacetimes, the lowest-order 
terms in the transformations do reduce to the well-known PM gauge transformations for a 
composite PM held. We also considered the bimetric model at the quadratic level and showed 
that it provides a unihed description of both linear PM theory as well as linearised conformal 
gravity, including in hat spacetime. We argued that the bimetric setup provides a natural 
and more promising framework for hnding a nonlinear PM theory, without a fundamental 
PM held, as compared to the usual Fierz-Pauli framework. 

Our analysis of the transformations as extensions of Weyl scaling is intimately connected 
to the g' -f-)- f interchange symmetry of the equations and breaks down in the massive gravity 
limit of the theory. Also the equations breakdown for the relevant bimetric theory in 3 
dimensions. Such two-derivative theories do not exist in other dimensions. Thus the analysis 
indicates the absence of such a symmetry in d = 3, as well as in massive gravity. 

If the 6-derivative transformation of the equations of motion found here can be extended 
to all orders, this would imply a local gauge invariance of the original equations. In the 
present approach going beyond the 6-derivative order is cumbersome, so the analysis is not 
conclusive. However, it is evident from our findings that the equations in this bimetric model 
possess a lot of hidden structure. A possibility is that a gauge symmetry exists and a better 
understanding of the structure of equations may enable one argue for its presence and even 
hnd its closed form expression. The model would then propagate six instead of seven degrees 
of freedom around any background and would become a ghost-free replacement for conformal 
gravity. In this case, imbedding the model in a larger setup should make the symmetry 
manifest. 

It is also possible that the equations do not have a gauge symmetry and our construction 
cannot be extended to all orders. Nonetheless, the results show that the bimetric framework 
is more powerful than the Fierz-Pauli framework in searching for a nonlinear PM theory. It is 
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also possible that a gauge symmetry exists in expected completions of the model. It is obvious 
that the bimetric model is the spin-2 analogue of the Proca theory in curved background. 
Hence, just like Proca theory, one may need additional degrees of freedom to obtain the 
model through a Higgs-like mechanism. This or other extensions may be better candidates 
for finding a new gauge symmetry. If such a nonlinear PM theory exists it would be very 
interesting to consider it as an alternative for conformal gravity in the proposal of [45, 46], 
which attempts to embed the Standard Model in a renormalisable theory of gravity. 

Finally, the symmetry we have found (at 6-derivative level) is an on-shell symmetry of the 
equations of motion, except for the quadratic theory, where it is also an off-shell symmetry 
of the action. In terms of the nonlinear fields, even the constant scalings of the metrics at 
the zero derivative level are not symmetries of the action, whereas they trivially keep the 
equations of motion invariant. Of course, off-shell symmetries of the action always imply 
invariance of the equations of motion, but the converse is not always true. 

5.2 Discussion of the counter-arguments 

Recently, some arguments have been made against the existence of a scalar gauge symmetry 
in the bimetric model with parameters (2.14). We have shown that this model goes beyond 
other similar constructions, in particular in displaying PM and conformal gravity behaviours 
in different phases, but we do not make any claims beyond what we have explicitly computed. 
It is also interesting to see that the model can evade the counter arguments presented so far 
and could suggest ways of improving them. First of all, since the transformations we construct 
are on-shell invariances, any off-shell analysis at the level of the action must also remain valid 
on-shell to be relevant. 


Massive gravity based arguments: The main argument for the absence of a gauge sym¬ 
metry in the bimetric model is based on the corresponding massive gravity limit. The massive 
gravity model, obtained in the a —)• oo limit [47] of the bimetric model (2.14) (for the unsealed 
variables g, /), was studied in [27-29]. It also admits proportional backgrounds g = c“^/ with 
arbitrary c^, where / is now a fixed Einstein-de Sitter metric with Kg = The lin¬ 

ear theory around this background is precisely the linear PM theory reviewed in section 3.1. 
Arguments were presented in [27-29] that, in this case, the linear PM symmetry does not 
extend to the nonlinear theory. It was also argued that since massive gravity can be obtained 
as a limit of bimetric theory, the absence of a PM symmetry in the massive gravity model 
implies the same for the bimetric model [27-29]. 

This argument assumes that the bimetric gauge transformations have an a ^ oo limit. 
The results in [20] and in the present paper indicate that this is not the case, as was explained 
in [31] in detail. We summarise the key points here. 

(i) Our derivation of the symmetry to sixth order heavily depended on the g' •(-)• f' 
interchange symmetry of the equations, which is destroyed in the massive gravity limit. 
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(ii) Our starting point was that the ^-equation of motion can always be rendered invariant 
by an appropriate transformation of /. In massive gravity, / is a fixed metric that cannot be 
varied (not even within a class of solutions of the Einstein equation for /). 

(iii) The expansion (4.10a) does not exist in massive gravity, since a term with 2n deriva¬ 
tives diverges as in the limit and there are infinitely many terms with increasing powers 

of a. Hence the massive gravity equations are not writable as the Bach equation plus correc¬ 
tions and the Weyl invariance at the lowest-order in the derivative expansion is not obtainable. 
Massive gravity therefore lacks the remarkable structure that exists in the bimetric equations. 

(iv) The gauge transformations Ag, Af obtained from (4.10) or (4.11) also involve an 
expansion powers of a and hence diverge in the massive gravity limit. The solutions that 
are obtained by applying these transformations also do not exit in the massive gravity limit. 
This indicates that these gauge transformations indeed disappear in the massive gravity limit. 
Of course, these statements are based on a perturbative expansion which, in principle, could 
sum up into a closed expression with a well defined a —)• oo limit. But, as pointed out in 
[31], solutions that are singular in the limit do exist. The key point is that massive gravity 
is a limit of bimetric theory around specific classes of solutions that need to be specified 
before taking the limit (a massive gravity limit of the entire bimetric theory does not exist). 
As discussed in [31], for any solution with a well-defined massive gravity limit, there exists 
another solution which becomes singular in this limit. 

If now a gauge transformation in nonlinear bimetric theory connects these two types of 
solutions, the perturbative expansion for the transformations will contain an infinite number 
of increasing powers of a (as is the case for the transformations constructed here), then clearly 
it cannot survive the massive gravity limit and the corresponding solutions are also absent 
after the limit has been taken. The form of the transformations derived here suggests that, 
generically, every gauge orbit contains one solution around which the bimetric theory has a 
massive gravity limit and the rest of the gauge orbit is invisible in the massive gravity limit. 

(v) So far we have argued that if the bimetric model (2.14) has a gauge symmetry of 

the type discussed here, the symmetry will not survive in the massive gravity limit. But 
this argument must exclude the PM symmetry of linearised massive gravity around EdS 
backgrounds. Indeed, in our analysis, the terms that diverge as a —>■ oo are for n >2. 

As mentioned earlier, on EdS backgrounds, Pf^uig) = both as well as 

their hrst variations vanish (4.12). Hence in this case the transformations do have a massive 
gravity limit, as should be the case. 

Symmetry algebra arguments: Some recent interesting works have studied the restric¬ 
tions imposed on nonlinear PM theories by the closure of the PM symmetry algebras [32, 35]. 
Reference [32] considered a general action for a PM field (j)^y interacting with a metric 
and analysed the theory to cubic order in At the lowest order, was assumed to 

be invariant and transformed as a linear PM field, but now in a general background. 
0/ii/-dependent corrections to the transformations were determined from the invariance of the 

^®An example is the flat bimetric solution in section 3.3 that is singular in the massive gravity limit 


- 25 ~ 




quadratic and cubic terms. On imposing the closure of the algebra it was found that 
whenever the PM field had a positive kinetic term, the cubic interactions were imaginary, 
while real cubic interactions were associated with a ghost PM field (as in conformal gravity). 
This lead to the conclusion that a unitary interacting theory with PM symmetry does not 
exist. [32] also emphasised that the action studied covered bimetric theory. If true, this 
would imply that the symmetry, at the 6-derivative level, found here is an artefact of the 
< 0 regime of bimetric theory and cannot be extended to higher orders. However, besides 
the analysis being off-shell, there is a potential caveat. To show that the PM interactions 
considered also covered bimetric theory, [32] used the relations, 

G ^1/ = fffj,!/ -|- ~ dflU f• (^- 1 ) 

and re-expressed the bimetric action (for rUg = nif) in terms of and Due to the 
symmetry S{g,f) = S{f,g) there will be no terms linear in (as assumed in [32]), but 
also there will be no cubic pg,y terms that were crucial in the analysis of [32]. Moreover, 
the identification of a PM field in bimetric theory is background dependent and the p^y in 
(5.1) can be identified with a PM field only around proportional backgrounds with (? = 1. 
Then, since a = 1, from the equations in section 4.4 it is evident that around proportional 
backgrounds, (corresponding to pgy = 0), the PM field does not transform. Note that in 
this case, the quadratic theory still has an accidental linear PM symmetry, but this does not 
follow from the nonlinear PM symmetry. The vanishing of the transformation for = 1 is 
completely consistent with the absence of the cubic p^y interactions. Hence, at least using 
the identifications (5.1) advertised in [32], the analysis there is avoided by the bimetric model. 

Another general study of PM algebra was carried out in [35]. The analysis involved 
only one dynamical spin-2 field but, nevertheless, seems related to our work. The authors 
considered general nonlinear extensions of the PM symmetry in dS and, by demanding closure 
of the transformation algebra, found a unique candidate transformation. The form of this 
symmetry, and the fact that it had a non-trivial flat space limit, led the authors to the 
conclusion that such an invariance could not be realised in any two-derivative single-field 
Lagrangian. These results are fully consistent with ours. In order to compare our setup to 
the single-field analysis of [35] we have to consider the single-field higher derivative equation 
(2.15). We found a candidate symmetry of those equations up to six-derivative order (in 
the transformation), which furthermore also had a flat space counterpart. It is obvious that 
the single-field equation (2.15) does not derive from any two-derivative single-field action. In 
fact, as a two-derivative theory they follow from the bimetric action containing both g^y and 
fgy whereas, when viewed as a single-field theory, they are higher-derivative equations which 
derive from a non-local action obtained after properly integrating out fgy. 

Finally, we should mention the constraint analysis of [33] in the d = 3 version of the 
present bimetric model given in [19]. The outcome, that the model in 3-dimensions does not 
posses PM symmetry is consistent with the fact that the expansions in (4.11) do not exist in 
the d = 3 theory. 
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It is interesting to see if a more refined version of the above analysis could make a more 


definitive statement about PM theory in the bimetric framework. 
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A Details of bimetric theory and its curvature expansion 
A.l Structure of action and equations of motion 

The elementary symmetric polynomials that enter the ghost-free bimetric action can be de¬ 
fined through the following recursion formula, 


rt — 1 




(A.l) 


k=0 


Note that 64 ( 5 ) = det(S') and 6^(5') = 0 for n > 4, if S' is a (4x4)-matrix. Explicit expressions 
can easily be written down. The contributions from the interaction potential to the bimetric 
equations of motions are given by [23], 




(A.2) 


in which we have used the definitions 


n 



(A.3) 


k=0 


Since the matrices S ^ and S with indices raised or lowered by either of the two metrics are 
symmetric, the same holds for the [y^"'^(S)]^i/. Note also the following identity. 


[Y^^\M)r, = 0 


(A.4) 


which holds for any 4x4 matrix M and follows directly from the Cayley-Hamilton theorem. 
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A.2 Obtaining the curvature expansions 

The explicit expression for th perturbative solution for obtained from the 5 ^ 1 / equa¬ 
tion (2.4a) reads [20], 


f/iu — 0, -\- 


2a? 

SIM' 


ZD I a^(si-l-2s2) 

if)? 


p2 

flL^ 


+ ^-^[le2iP)g^u-PP^.u]+0 


p3 


(A.5) 


where g? = m^jm^ and all indices on the right-hand side are contracted with g^y. The Sn 
are particular combinations of bimetric parameters given by, 



(A. 6 ) 


The value of the coefficient a? is obtained by solving the equation sq = 0 for a (where we 
only consider solutions that result in real values for a^). The solution for g^j^^ obtained from 
the /^jy equation (2.4b) has the very similar form, 




si/i- 


z3 I a^(si-|-2s2) 132 


fiu 


I 2o?S2 

PW 


ke2{P)U. - PP^ 


flU 


+ 0 


p3\ 

W 


(A.7) 


where P^u = Pfj,u{f) is the Schouten tensor for Indices are now contracted with and. 


= 


4 

m 

a^rrig 



(A. 8 ) 


with a being a solution to the polynomial equation §0 = 0 . 

We can use the expression (A.5) for to eliminate it from its own equation. For 
instance, inserting (A.5) into the Einstein tensor for results in, 


GAf) = GAd) - ^ - VPV^Pp, - VPV.Ppp 


+ WPp, - gp, [P^^Pa0 + \P^ 
The contributions from the interaction potential read as, 
f7 — n 4- 1 c 4- 20 p 

c? a^a^si 


+ 0 


p3 


) ■ (A.9) 




P 'a?o?sffj? 


XlPPPp, + X2PPpu + \gAx^P"^PaP - X2P^)\ + O (^) , 

in which the coefficients Xn are given by, 

xi = 2s^ -|- f4(si -|- 2 S 2 ), X2 = —3s^ — 2 S 2 GI , 3:3 = 3s^ — 252^7 , 


(A.IO) 


(All) 


and we have dehned, 


n = a/3i -|- 2)0? (52 + . 


(A. 12 ) 
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Combining the above results, the entire equation of motion becomes the following higher 
derivative equation for 






-'flu 


I 

£j2Q,2^a^4 


(si + 2S2)P/Pp. - 2s2PPp. - f 5/4. - P^) 




3PPp, - 2P/Pp, - -5m. (P2 - P“^P„;3) 


+ 0 


p3 

7^ 


= 0 . 


(A. 13 ) 


Here, some of the four-derivative terms have been collected into the Bach tensor [1], 

Bfiu = ~'^'^Pfiu — VpVj/P^ + VpVpP^ + 'S/p's/— 2Pp^Ppy + ^gpuP^'^Ppa ; (A.14) 
which is invariant under local Weyl transformations of the metric. 


B Four-derivative terms in gauge transformations 

Here we outline the computation of the four-derivative terms generated by lower orders in 
the symmetry transformation. 

B.l Higher-curvature equations 

We start by presenting the explicit expressions for the higher-curvature expansions of the 
equations of motion in the PM bimetric model up to 10th order in derivatives. The equations 
read, 


OO 




II 

1 

+ ^pPp,u 

n=2 



(B.la) 
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(B.lb) 

Here, the 

higher-derivative functions 

7 p?r^ for n < 5 are of 

the form,^^ 
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(B.2) 

^( 6 ) 
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= (136? - § 62 ) P,l 

- (I4e? 

- 76162 - 663 ) PI^ 







-h {e\ -h 86?62 + §6? - 

- 136163 

- 864 ) Pm. 



- (5efe2 - |eie| + bejes + 56263 - 226164) ■ (B. 5 ) 

have made use of the Cayley-Hamilton theorem (A.4) to reduce the tensor powers in 7 ® and 7 ^ii°^. 
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In this expression all of the are elementary symmetric polynomials of the Schouten tensor 

P't = 9^PPp.{9). 

For future reference we also present the following linear variations of the four-derivative 
terms, 



2P/P/ - P,,PP^ - \e2{P)5P^5l + \9puPP^^ - \ 9 puP\P 


Act 


^ 9 , 


pa 
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1 
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Ppu9'^^ + P5';.K - 2P^P5Z - 2P/6; - -9puP9'^^ + ^9puP'^^ 


6 P 


per ? 


(B.6) 


where the linearised Schouten tensor reads as, 


5Pf,ty = P^y[g 5g\ - Ppv[9] 

= -\ [VpVuSg - - VPVJgpp + V‘^5gp,] + ^gp, [V‘^5g - VPV^dg^p] 


2 

+ ^9pu5gP^Ppa - ^P^p 


^9py Q9py^9 


(B.7) 


and we have. 


6en{P) = ^i-l)'^en-kiP) \[P^rSgp^ - [P^-^rSPp^ 


k=l 


(B.8) 


These will be useful for computing the contribution to the four-derivative terms in the gauge 
transformations. 


B.2 Vanishing of the four-derivative contributions 

Consider arbitrary variations expanded in powers of derivatives up to 4**^ order, 

‘^g'py = ‘^{0)9'py + ^(2)9'pu + ^(4,)9'pu + ■ ■ • , 

^f'pu = ‘^{0)fpu + ^{2)f'pu + ^{i)fp.y + • • • • (B.9) 

The terms with zero and two derivatives were determined earlier, 

^{0)5^^ = (1^9'pu > ^i2)9pu = -‘^p(I^Ppy - 2pVpdu(t) , 

^(0)f'pu = -(l^g'pu , ^{2)fpu = ‘2^p(t>Ppu - 2pVpdy(l). (B.IO) 

Using (B.lb) to replace = —f'^y + ^pPpv + • • •, the transformations of f'^y can also be 
expressed in terms of f'py, 

^{0)fpy = (t>f'py , ^{2)fpu = -‘2^P(t>Ppu “ 2pVpdyCj) . (B.ll) 
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Demanding the invariance of (B.la) at the four-derivative level and using the results from 
appendix B.l, we then obtain the following relation among the transformations, 
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- 2P/VpV,</. - 2P/VpVpcP - kg'PV^^ + P^'^VpV 


(B.12) 


where all quantities on the right hand side are given in terms of g'^^. We now use (B.lb) to 
express the right hand side of (B.12) as a functional of f'^^, 

^{0)f'pu + ^{2)f'pu + ^(4)/pi/ 
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(B.13) 


where now, on the right-hand side, all quantities (curvatures, derivatives, index raising, etc.) 
are understood to be defined with respect to and in the transformations g'^^^ is a function 
of /pi/, 


^{2)9pv ^{2)9pu\gi=_fij^4^pp ) 

Furthermore, we have introduced 
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(B.14) 


(B.15) 


p-* 1/ W V i// ^ V / pj/ 

Next, we try to impose the interchange symmetry of the transformations. According to our 
discussion in section 4.3, this is only possible if the contributions to the transformations that 
are generated from lower orders vanish. Indeed, if we demand, 


^(4) /pi/ ^{A)9pi/\gl^p , 

after a lengthy but straightforward calculation we find that this condition becomes, 
A(4)5;, = A(4)5;, + 4p2 [VpVpVi/S^c/- + VPV.Vpdpcj) - VpVi/V^,^ - V^Vp^i/c^ 


(B.16) 
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Clearly, this leaves undetermined and instead provides a consistency check on the 

possible existence of a symmetry to this order in derivatives. Obviously, the condition is 
satisfied for constant gauge parameter, but showing that the extra terms vanish for a general 
function cj) requires a bit more work. 

Some identities that are useful for our purposes at this point are,^® 

VpV^V,V^</> = 

+ P/V,Vp</. + + VpP.pVPc /), (B.18) 

V^V.VpVpc/. = + VPR^p^^V^cj ), (B.19) 

VpVaVPV^cf) = + PP^VpVaCp + \P^^(P + ^VpPV^c/). (B.20) 

Using these identities, we can write the condition (B.17) as, 

VpR^,/V^^ + VpP.^V^cj) - V„Pp,V-ct> + = 0 . (B.21) 

Now, in order to show that this is always true, consider the (contracted) second Bianchi 
identity. 


VpPp,/+ VpP,, - = 0 , (B.22) 

or equivalently, using Rp^ = Ppy + ^QpuP, 

P^pau^ + ^pPyo — ^aPpu + -Quct'^ pP — -Qpv^aP = 0 . (B.23) 

Using this, it is obvious that (B.21) and hence (B.17) are identically satisfied for any func¬ 
tion cj). This proves that there are no four-derivative terms generated from lower orders in the 
gauge transformations and hence the condition for the existence of an invariance is satisfied 
at the first nontrivial order. 
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